Towards the observation of Hawking radiation in Bose— Einstein condensates. 
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Acoustic analogues of black holes (dumb holes) are generated when a supersonic fluid flow entrains 
sound waves and forms a trapped region from which sound cannot escape. The surface of no return, 
the acoustic horizon, is qualitatively very similar to the event horizon of a general relativity black 
hole. In particular Hawking radiation (a thermal bath of phonons with temperature proportional 
to the "surface gravity") is expected to occur. In this note we consider quasi-one-dimensional su- 
personic flow of a Bose-Einstein condensate (BEG) in a Laval nozzle (converging-diverging nozzle) , 
with a view to finding which experimental settings could magnify this effect and provide an observ- 
able signal. We identify an experimentally plausible configuration with a Hawking temperature of 
order 70 n K; to be contrasted with a condensation temperature of the order of 90 n K. 
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Acoustic analogues of black holes are formed by super- 
sonic fluid flow [0, H] . The flow entrains sound waves and 
forms a trapped region from which sound cannot escape. 
The surface of no return, the acoustic horizon, is qual- 
itatively very similar to the event horizon of a general 
relativity black hole; in particular Hawking radiation (in 
this case a thermal bath of phonons with temperature 
proportional to the "surface gravity" ) is expected to oc- 
cur ||l|, [2). There are at least three physical situations 
in which acoustic horizons are known to occur: Bondi- 
Hoyle accretion , the Parker wind 0| (coronal outflow 
from a star), and supersonic wind tunnels. Recent im- 
provements in the creation and control of Bose-Einstein 
condensates (see e.g., [^) have lead to a growing interest 
in these systems as experimental realizations of acoustic 
analogs of event horizons. In this note we consider super- 
sonic flow of a BEC through a Laval nozzle (converging- 
diverging nozzle) in a quasi-one-dimensional approxima- 
tion. We show that this geometry allows the existence of 
a fluid flow with acoustic horizons without requiring any 
special external potential, and we study this flow with 
a view to finding situations in which the Hawking effect 
is large. We present simple physical estimates for the 
"surface gravity" and Hawking temperature . 

While writing up this paper, we have found that fiuid 
flow in a Laval nozzle geometry has also been considered 
in reference . The occurrence of the same Laval nozzle 
geometry is the only significant point of overlap of our 
work with j?] , as that paper only deals with a classical ef- 
fect, related to the Hawking effect, but does not consider 
the quantum physics of Hawking radiation itself. (See ||] 
for a previous discussion along these lines.) 



Laval nozzle: A general problem with the realization 
of acoustic horizons is that most of the background fluid 
flows so far studied seem to require very special fine- 
tuned forms for the external potential in order to be re- 
alized (see e.g., the Schwarzschild-like geometry in refer- 
ence Q). In this respect a possible improvement toward 
the realizability of acoustic horizons is the construction 
of a trap which "geometrically constrains" the flow in 
such a way as to replace the need for a special exter- 
nal potential. An example of such a geometry is the so 
called Laval nozzle. In particular we shall consider a pair 
of Laval nozzles; this provides a system which includes a 
region of supersonic flow bounded between two subsonic 
regions. 




FIG. 1: A pair of Laval nozzles: The second constriction is 
used to bring the Quid How back to subsonic velocities. 

Consider such a nozzle pointing along the z axis. Let 
the cross sectional area be denoted A{z). We apply, with 
appropriate modifications and simplifications, the calcu- 
lation of |9| . The crucial approximation is that transverse 
velocities (in the x and y directions) are small with re- 
spect to velocity along the z axis. Then, assuming steady 
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flow, we can write the continuity equation in the form 

p{z) A{z) v{z) = J; J — constant. (1) 

The Euler equation (including for the moment possible 
external body forces d<f>/dz, and internal viscous friction 
fv) reduces to 



dv dp d<i> 

dz dz dz 
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Finally, we assume a barotropic equation of state p = 
p{p)^ and define X' = dX/dz. Then continuity implies 
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Defining the speed of sound by = dp/ dp, and elimi- 
nating p' between these two equations yields a form of 
the well-known "nozzle equation" 
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The presence of the factor — v"^ in the denominator is 
crucial and leads to several interesting physical effects. 
For instance, if the physical acceleration is to be finite at 
the acoustic horizon, we need 
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This is a fine-tuning condition relating the shape of the 
nozzle (the amount of convergence/divergence at the 
horizon) to the external body force and specific friction. 
Experience with wind tunnels has shown that the flow 
will attempt to self-adjust (in particular, the location of 
the acoustic horizon will self-adjust) so as to satisfy this 
fine tuning. We could now calculate the acceleration of 
the fiuid at the acoustic horizon by adopting L'Hopital's 
rule. However, it is more useful to consider the "surface 
gravity" defined by the limit of the quantity S 
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It is this combination g, rather than the physical accel- 
eration of the fluid a, that more closely tracks the notion 
of "surface gravity" , and it is the limit of this quantity 
as one approaches the acoustic horizon that enters into 
the Hawking radiation calculation . To calculate the 
limit gH we need to use 
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This implies that the flne-tuning (g) used to keep a finite 
at the acoustic horizon will also keep g finite there. In 
particular, we can use L'Hopital's rule at the acoustic 
horizon to evaluate 
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This result, which we quote for its generality, is actually 
more complicated than we really need. 

Free flow: The Laval nozzle corresponds to the special 
case in which there is no external body force = 0. 
Additionally, when considering a superfluid flow, there is 
no viscous friction = 0. Then the nozzle equation (|^) 
reduces to 
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Regularity now requires the much simpler fine-tuning 
condition that A' — &t the horizon. That is, the acous- 
tic horizon occurs at a point of minimum area; exactly 
the behaviour which is physically seen in a Laval noz- 
zle. (That a horizon cannot form at a maximum of the 
cross sectional area is established below.) Now apply the 
L'Hopital rule at the horizon (using the fact that at this 
point A' = 0) 
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But from the specialization of (|^) we now see 
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Very similar formulae hold for g^'- 
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The first factor is of order cj^/R, with R the minimum 
radius of the nozzle, while the second and third factors 
are square roots of dimensionless numbers. This is in 
accord with our intuition based on dimensional analy- 
sis [|[ ^. If A" < 0, corresponding to a maximum of the 
cross section, then gn is imaginary which means no event 
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horizon can form there. The two signs ± correspond to 
either speeding up and slowing down as you cross the 
horizon, both of these must occur at a minimum of the 
cross sectional area A" > 0. (If the flow accelerates at 
the horizon this is a black hole horizon [future horizon] ; if 
the flow decelerates there it is a white hole horizon [past 
horizon]. See Figure 1.) If the nozzle has a circular cross 
section, then the quantity A'l^ is related to the longitu- 
dinal radius of curvature Rc at the throat of the nozzle, 
in fact A'^ = nR/Rc. 

Bose-Einstein condensate: The use of BECs as a work- 
ing fluid for acoustic black holes has been advocated by 
Garay et al. |ll[ (see also [Q for a discussion of plausible 
extensions to that model). The present note can be in- 
terpreted as a somewhat different approach to the same 
physical problem, side-stepping the technical complica- 
tions of the Bogoliubov equations in favour of a more 
fluid dynamical point of view. For a standard BEG 
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So we have, rather simply 
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This implies, at a black hole horizon [future horizon], a 
Hawking temperature [Q, 0, nfl 
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Ignoring for now the issue of gray-body factors (they are 
a refinement on the Hawking effect, not really an essential 
part of the physics) , the phonon spectrum peaks at 
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This extremely simple result relates the Hawking emis- 
sion to the physical size of the constriction and a factor 
depending on the flare-out at the narrowest point. Note 
that you cannot permit A'^^ to become large, since then 
you would violate the quasi-one-dimensional approxima- 
tion for the fluid flow that we have been using in this note. 
(There is of course nothing physically wrong with vio- 
lating the quasi-one-dimensional approximation, it just 
means the analysis becomes more complicated. In par- 
ticular, if there is no external body force and the viscous 
forces are zero then by slightly adapting the analysis of |^ 
the acoustic horizon [more precisely the ergo-surface] is 
a minimal surface of zero extrinsic curvature.) The pre- 
ceding argument suggests strongly that the best we can 
realistically hope for is that the spectrum peaks at wave- 
length 
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You can (in principle) try to adjust the equation of state 
to make the second factor in ( [l5| ) larger, but this is un- 
likely to be technologically feasible. 

(Note that this is the analog, in the context of acoustic 
black holes, of the fact that the Hawking flux from gen- 
eral relativity black holes is expected to peak at wave- 
lengths near the physical diameter of the black hole, its 
Schwarzschild radius — up to numerical factors depend- 
ing on charge and angular momentum.) 

Discussion: It is the fact that the peak wavelength of 
the Hawking radiation is of order the physical dimensions 
of the system under consideration that makes the effect so 
difficult to detect. In particular, in BEGs it is common 
to have a sound speed of order 6 mm/s. If one then 
chooses a nozzle of diameter about 1 micron, and a flare- 
out of A'Ij Ri 1, then Th ~ 7 n K. Gompare this to the 
condensation temperature required to form the BEG 



^condcns 



90 nK. 
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We see that in this situation the Hawking effect, although 
tiny, is at least comparable in magnitude to other relevant 
temperature scales. Moreover recent experiments indi- 
cate that it is likely that these figures can be improved. 
In particular, the scattering length for the condensate can 
be tuned by making use of the so called Feshbach reso- 
nance jl^ . This effect can be used to increment the scat- 
tering length; factors of up to 100 have been experimen- 
tally obtained [Q. Therefore the acoustic propagation 
speed, which scales as the square root of the scattering 
length, could thereby be enhanced by a factor up to 10. 
This suggests that it might be experimentally possible to 
achieve c// ~ 6 cm/s, and so 



w 70 n K, 



(26) 



which places us much closer to the condensation temper- 
ature. The speed of sound can also be enhanced by in- 
creasing the density of the condensate (propagation speed 
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scales as the square root of the density) . In all of these 
situations there is a trade-off: For fixed nozzle geometry 
the Hawking temperature scales as the speed of sound, 
so larger sound speed gives a bigger effect but conversely 
makes it more difficult to set up the supersonic flow. 

So if the Hawking effect can be experimentally realized 
in these situations, it may be sufficiently large to disrupt 
the condensate configuration, or even the condensate it- 
self, providing in this way a clear signal. Let us elaborate 
this point: The Garay et al. analysis ||ri| shows that it 
should be possible to create classically stable BEG con- 
figurations with the presence of acoustic horizons. In 
their particular analysis they found that these classically 
stable configurations are surrounded by unstable regions; 
(at this stage, we don't know how general that result is). 
Once the system is engineered to be in a classically stable 
configuration one can look for the purely quantum effect 
of Hawking emission. The power loss due to Hawking 
radifiLluu would be. 



lated: via the speed of sound, the nozzle radius, the equa- 
tion of state, and the degree of flare-out at the throat. In 
this manner, we have identifled an experimentally plausi- 
ble configuration with a Hawking temperature of order 70 
n K; to be contrasted with a condensation temperature 
of the order of 90 n K. 
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Numerically (including the effect of the Feshbach res- 



ona wcc), — the emitted power is extremely small P — 
IQ— w but it should be noted that some sizable frac - 
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tion of the Hawking phonons will be sufficiently energetic 
to knock atoms out of the condensate phase. 

Additionally, we want to stress that the present anal- 
ysis of the BEG is purely "hydrodynamic" (superfiuid 
approximation of the BEG), and does not seek to deal 
with the "quantum potential" (see e.g., |l2|). The lat- 
ter is responsible for the fact that the dispersion relation 
for perturbations in the BEG is modified at high mo- 
menta in such a way as to recover "infinite" propagation 
speed (this is the so-called Bogoliubov dispersion rela- 
tion Ijl^, |l^). This issue has relevance to the so-called 
trans-Planckian problem (which in this BEG condensate 
context becomes a trans-Bohrian problem). Fortunately 
it is known, thanks to model calculations in field theo- 
ries with explicit high-momentum cutoffs, that the low 
energy physics of the emitted radiation is largely insen- 
sitive to the nature and specific features of the cutoff. 
However, because of this high-frequency "superluminal" 
dispersion, one can have additional hopes of detecting 
a signal because Hawking radiation could provide a new 
type of instability, disrupting the classically stable config- 
uration by leading to a run-away production of phonons 
(see the discussion of "black hole lasers" in jl^). 

To summarize what we have learned: The present note 
complements the analysis by Garay et al. |l^, in that it 
provides a rationale for simple physical estimates of the 
Hawking radiation temperature without having to solve 
the full Bogoliubov equations. Additionally, the current 
analysis provides simple numerical estimates of the size 
of the effect and identifies several specific physical mecha- 
nisms by which the Hawking temperature can be manipu- 
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